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Description 

BACKGROUND OF THE INVENTION 
5 Technical Field 

The present invention is directed to filtering well logging data. More particularly, the present invention is 
directed to filtering out unwanted components contained in acquired sonic well logging data, e.g., reflected 
wave components, converted wave components, casing arrivals, tool arrivals and/or noise. 
w The present invention is further directed to estimating a reflection coefficient of a bed boundary in the 

formation, identifying the lithology of the formation traversing the borehole, and compressing the sonic data, 
based on the acquired sonic data. 



Background Information 

15 

Turning now to Figure 1 , a schematic diagram of a logging operation is shown. Tool, or sonde, 10 for ac- 
quiring sonic data is located in borehole 11 penetrating earth formation 12. Optionally, the borehole wall may 
have casing 13 cemented thereto, e.g., in a production well. The sonde is preferably lowered in the borehole 
by armored multi-conductor cable 14 and slowly raised by surface equipment 15 over sheave wheel 16 while 
20 sonic data measurements are recorded. The depth of the tool is measured by depth gauge 17, which measures 
cable displacement. 

. Sonde 1 0 acquires sonic data by emitting an acoustic pulse and recording its return waveform. The sonde 
comprises at least one sonic source, or transmitter, and at least one detector, or receiver. The source typically 
produces a pulse upon excitation. The pulse travels through the casing and formation, and back to the sonde 
25 where U is detected by the receiver(s) thereon. The return waveforms can be analyzed by the sonde in situ, 
analyzed by data processor 18 at the surface, or stored, either in the sonde or at the site, for analysis at a 
remote location. In the preferred embodiment, the return waveform data is transferred to data processor 18 
by cable 14 for analysis at a remote location. 

Sonic data acquired in this manner is typically displayed on a chart, or log, of waveform amplitude over 
30 time versus depth. With reference to Figure 2, a representation of sonic data in log format is shown. The data 
was recorded at 200 successive depth locations, the depth interval consisting of various formations, ranging 
from very hard limestone to very soft shaley sand. 

Ideally, the waveform data records the arrival of compressional (P) waves, shear (S) waves and Stoneley 
(ST) waves. However, the P, S and/or ST wave components often contain undesirable components which de- 
35 grade and/or mask the desired wave components. Examples of undesirable components include reflect d 
and/or converted waves, casing arrivals, tool arrivals and noise. 

The P, S and/or ST wave components are often contaminated by reflected and/or converted waves when 
the receiver is close to a bed boundary. Reflected waves are waves which are emitted by the source as P, S 
or ST waves, reflected by bed boundaries, and received as P, S or ST waves, respectively. Converted waves, 
40 on the other hand, are waves which are emitted by the source as P, S or ST waves, reflected by bed boundaries, 
and received as something other than P, S or ST waves, respectively. For example, a P wave which is "con- 
verted" to an S wave, and so recorded, is a converted wave. 

Several examples of reflected and converted waves are marked at A through I in Figure 2. As can be seen 
in Figure 2, the reflected and converted waves occur at the same time as the desired waveform arrivals. Ad- 
45 ditionally, the frequency content of the reflected and converted waves overlaps that of the P, S and/or ST wave 
components. 

Aside from P, S and ST wave component arrivals from the formation, the acoustic pulse travels through 
the steel casing, as well as the tool itself. The waveforms therefrom are commonly referred to as casing and 
tool arrivals, respectively. The point in time at which a waveform component arrives and is subsequently de- 
50 tected by the sonde's receiver(s) is commonly referred to as the waveform component's arrival time. 

A steel casing is typically cemented to the borehole wall in production wells. As is well known, the cemented 
casing isolates the various water, gas and oil bearing zones from each other, thereby maintaining zone integrity. 
When the cement bond is good, the casing is substantially transparent to the sonic tool. Thus, the sonic tool 
is capable of logging through the casing and acquiring the formation wave components. However, when the 
55 cement bond is deficient or deteriorated, a strong wave travels through the casing, subsequently detected by 
the receiver(s). The casing arrival can mask the formation arrivals which have a smaller amplitude. 

The tool itself also provides a path from transmitter to receiver. The tool's arrival time can appear at the 
same time as the formation arrival, based on tool design. Further, the amplitude of the tool arrival can be several 
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times that of the formation arrival. As a result, the formation wave is masked by the tool arrival. As appreciated 
by those skilled in the art, the formation wave is desired for many reasons, eg, to calculate the formation slow- 
ness. 

Additionally, the recorded waveforms are likely to be contaminated by noise, eg, geological, external, elec- 

5 tronic and/or quantization. Geological noise is generated by the formation, eg, from fracture dev lopment. Ex- 
ternal noice is due to interference from external sources, eg, road traffic, rig activity and the like. Electronic 
noise is produced by the electronic components of the sonde, eg, due to component reaction to thermal fluc- 
tuations, cross-talk or shot-noise. Quantization noise is the result of waveform degradation inherent in digitising 
an analog signal, such as the waveforms acquired by the sonic detectors. Since noise tends to degrade and 

10 mask the acquired data, noise is undesirable. 

Numerous filtering techniques are known for removing undesired components from waveforms. For ex- 
ample, the waveform can be filtered in either the temporal or the spectral domain. Such techniques, how ver, 
cannot be appl ied to sonic data to remove undesired components caused by reflected and/or converted waves, 
casing arrivals, tool arrivals and/or noise. In order for temporal or spectral filtering to be effective, there must 

15 exist a time or frequency separation, respectively, between the desired and undesired components. Such sep- 
aration does not exist between the desired and undesired components in the acquired waveform. 

Additionally, the waveforms can be filtered using velocity filtering techniques, such as that described for 
sonic well logging waveforms in the paper entitled "Velocity Filtering of Acoustic Well Logging Waveforms" by 
Kai Hsu and Thomas L Marzetta, published in the IEEE TRANSACTIONS ON ACOUNSTICS, SPEECH AND 

20 SIGNAL PROCESSING, vol. 37, no. 2, February 1989, pages 265 and 273. Reflected waveforms appear at a 
receiver after the direct waveform, typically having similar frequency characteristics, albeit different velocities 
(ie, arrival times). In order to separate the two, a velocity filter is employed which separates arrivals based on 
their different velocities. 

However, velocity filtering does not filter out noise components. Additionally, conventional filtering tech- 
25 niques do not allow the identification of lithology, let alone allow for data compression. 

Another filtering technique, for suppressing noise in seismic data, is decribed in the paper "Coherent Noise 
Suppression in Seismic Data using Eigenvalue/Eigenvector Decomposition", co-authored by William J Done 
and R Lynn Kirlin and published in the Proceedings of the 1988 International Conference on Acoustics. Speech 
and Signal Procesing (ICASSP 88), pages 908 et seq. 

30 

SUMMARY OF THE INVENTION 

Accordingly, it is an object of the present invention to substantially filter out undesired components, eg, 
reflected wave components, converted wave components, noise, casing arrivals, and/or tool arrivals from ac- 
35 quired sonic data. 

It is also an object of the present invention to estimate a reflection coefficient of a bed boundary in the 
formation for at least one wave component of the acquired sonic data. 

It is a further object of the present invention to identify the lithology of the formation traversing the borehol 
from which the sonic data was acquired. 
40 It is yet a further object of the present invention to compress the acquired sonic data. 

Four different aspects of the invention are defined hereinafter in claims 1 , 5, 8 and 1 2 respectively. 

Sonic data is recorded downhole, typically in digital format, by any conventional technique well known in 
the art. The acquired data set comprises a plurality of waveforms, each having compressional (P), shear (S) 
and/or Stoneley (ST) wave components. These componnts are operated upon individually by the method of 
45 the present invention. In the preferred embodiment, the P, ST and then the S wave components are filtered. 

In order to remove reflected wave components, converted wave components and/or noise, the position 
of the P wave component's arrival time for each waveform in the data set is detected. Based on this information, 
the P wave component is isolated from the rest of the waveform by a time window, predetermined to capture 
the entire P wave component. 
50 In the first embodiment, the n digitized samples of each P wave component are treated as an n-d imentional 

vector. These vectors can be plotted mathematically in an n-dimensional space. Given the similar nature of P 
wave components, the vectors will have a substantially similar primary direction. Fluctuation about and pro- 
jection on this primary direction is a function of the formation composition traversed by the borehole, as well 
as undesired components recorded in the wav forms. 
55 Karhunen-Loeve rotates the coordinate system of the n-dimensional space to a new coordinate system 

whos directions, characterized as eigenvectors, are aligned with the primary directions of the wav compo- 
nent vectors. The eigenvectors are ordered such that the first eig nvector represents the wave component's 
main direction, the second eig nvector represents the wave component's next main dir ction, et cetera. By 
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projecting the waveforms onto the principal eigenvectors, the reflected wave components, converted wave 
components and/or noise component of each waveform is substantially eliminated. 

In the second embodiment, an alternative method is employed for removing reflected wave components, 
converted wave components and/or noise. The n samples of each P wave component are averaged to obtain 
5 a template of the P wave components. The template is normalized, and this normalized template closely ap- 
proximates the first eigenvector discussed above. An approximation to the second eigenvector is preferably 
obtained by taking the Hilbert transform of the normalized template. 

Theoretically, there are n eigenvectors, corresponding to the original n-dimensional space. Practically, 
however, given the similar nature of the P, S and ST wave components, fluctuation about the first eigenvector 
10 is substantially limited. I have found that only the first few eigenvectors are needed to substantially describe 
the position of the wave component vectors. Thus, using only the first eigenvector, or an approximation thereof, 
about 86% of the wave components are captured. Using the first and second eigenvectors, or approximations 
thereof, about 92% of the wave components are captured; using the first three, about 96%. 

In both embodiments, the first eigenvector, or its approximation, is correlated to the individual wave com- 
15 ponent, e.g., P, for each waveform, thereby obtaining a correlation factor. Multiplying the correlation factor by 
the first eigenvector, or its approximation, yields a correlated eigenvector. This correlated eigenvector is sub- 
tracted from the wave component. The residual wave component includes undesired components, e.g., reflect- 
ed wave components, converted wave components, and/or noise, as well as residual wave components. In 
order to reduce the amount of wave component in the residual waveform, the second eigenvector, or its ap- 
20 proximation, is correlated to the residual waveform and subtracted therefrom. This correlatio/subtraction proc- 
ess can continue for all subsequent eigenvectors, based on desired resolution. 

The residual waveform contains the undesired components, such as reflected waves, converted waves, 
and/or noise. This residual P wave component waveform is subtracted from the acquired (original) wave com- 
ponent waveform, thereby producing wave components substantially free of these undesired components. 
25 This method, repeated for all wave components, thereby produces better slowness logs. 

In order to remove the undesired components of tool arrivals and casing arrivals, the tool arrival (casing 
arrival) is detected, and its first eigenvector, or its approximation, is determined. A correlation factor is calcu- 
lated and multiplied by the first eigenvector, or its approximation, yielding a correlated first eigenvector. The 
correlated first eigenvector is subtracted from the waveform. Optionally, these steps may be repeated a plur- 
30 ality of times to account for the second through nth eigenvectors. Due to the consistent nature of the tool arrival 
(casing arrival), once subtraction has occurred, the remaining waveform is substantially free of tool arrivals 
(casing arrivals). In the preferred embodiment, two eigenvectors, or their equivalent, are employed for removing 
casing arrivals, and one eigenvector, or its equivalent, is preferably employed for removing tool arrivals. 

In the preferred embodiment, the tool and casing arrivals are removed prior to removal of the other unde- 
35 sired components from the acquired waveforms. 

In another aspect of the present invention, the reflection coefficient of a bed boundary in the formation 
can be estimated for each wave component, once the wave components are substantially free of the above- 
mentioned undesired components. The energy content of the reflected wave component can be obtained from 
the residual waveform. The energy content of the transmitted wave component can be obtained from the f il- 
40 tered wave component. The reflection coefficient is estimated by dividing the energy content of the ref I cted 
wave component by the energy content of the transmitted wave component. 

In a further aspect of the present invention, the lithology of the formation traversing the borehole can also 
be identified. The P wave component, for example, can be expressed in a two-dimensional Karhunen-Loeve 
plane. The plane's abscissa (x-axis) is preferably the projection of the respective P wave component on the 
45 first eigenvector (obtained from the above-calculated factor correlating the P wave component vector to the 
first eigenvector or its approximation). The plane's ordinate (y-axis) is preferably the projection of the respec- 
tive P wave component on the second eigenvector (obtained from the above-calculated factor correlating the 
P wave component vector to the second eigenvector or its approximation). 

These points are then plotted in the Karhunen-Loeve plane. I have found that these points tend to cluster 
so based on lithology type. For example, points obtained from wave components traveling through limestone tend 
to congregate in one area of the plane, while points obtained from wave components traveling through shaley 
sand tend to congregate in another area. 

Given a plurality of wave compon nts depicted in this manner, a database can be established for identi- 
fying the lithology of the formation traversing the borehole given its coordinates in the Karhunen-Loeve plane. 
55 Given the database and the waveform's eigenvectors, the lithology of the formation through which the wave 
traveled can be identified. 

In yet a further aspect of the present invention, the data set can also be compressed, employing the ei- 
genvectors. Recalling that Karhunen-Loeve rotates the n-dimensional coordinate system to a new coordinat 
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system, the wave component vectors can be expressed in terms of the primary eigenvectors. For example, if 
the 4th through 50th eigenvalues are substantially equal to zero, each wave component vector can substantially 
be described using only their projection on the first three eigenvectors. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 illustrates a schematic diagram of a logging operation. 
Figure 2 illustrates acquired sonic logging data in log format. 

Figures 3a and 3b illustrate a first embodiment of the present invention for substantially removing reflected 
10 wave components, converted wave components and/or noise from acquired sonic data employing a Karhunen- 
Loeve transformation. 

Figure 4 shows the sonic data of Figure 2 such that their respective P wave component arrival times are 
aligned. 

Figure 5 shows the sonic data of Figure 2 wherein the reflected wave components, converted wave com- 
15 ponents and/or noise, calculated from the first eigenvector, are removed therefrom. 

Figure 6 shows the sonic data of Figure 2 wherein the reflected wave components, converted wave com- 
ponents and/or noise, calculated from the first and second eigenvectors, are removed therefrom. 

Figure 7 illustrates a second embodiment of the present invention for substantially removing reflected 
wave components, converted wave components and/or noise from acquired sonic data employing an approx- 
20 imation of the Karhunen-Loeve transformation of Figure 3. 

Figure 8 illustrate a first embodiment of the present invention for substantially removing casing and tool 
arrivals from acquired sonic data employing a Karhunen-Loeve transformation. 
Figure 9 illustrates an acquired sonic data set having multiple tool arrivals. 
Figure 10 illustrates the data set of Figure 9 with the tool arrivals removed therefrom. 
25 Figure 11 illustrates the identification of formation lithology based on a plane formed from the projections 

of the P wave components on the first and second eigenvectors. 

DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 

30 Turning now to Figures 3a and 3b, the first embodiment of the method of the present invention for sub- 

stantially removing reflected wave components, converted wave components and/or noise from acquired sonic 
data employing a Karhunen-Loeve transformation is described. 

At step 301 , the sonic data is acquired. In the preferred embodiment, the waveforms are digitized using a 
sampling rate of 10 us. Other sampling rates, depending upon desired accuracy and the Nyquest criterion, will 

35 be obvious to one skilled in the art. Methods of acquiring sonic data are well known to those skilled in the art, 
for example, as shown with reference to Figure 1 . Sonic data is typically acquired from substantially the entire 
length of the borehole. The data set shown in Figure 2 comprises 200 waveforms, each having P, S and ST 
wave components. 

Preferably, a 500 microsecond (us) time window is employed with a sampling rate of 10 us, thereby ob- 
40 taining 50 sampled points for the P wave components of each waveform. 

It is to be understood, however, that the dimensions of the data set used herein is for illustrative purposes 
only. Additionally, although the data set used herein includes P, S and ST wave components, it is to be under- 
stood that a sonic data set having any number of wave components can be used in conjunction with the present 
invention. 

45 As described hereinbelow, the method of the present invention is applied to remove reflected wave com- 

ponents, converted wave components and/or noise, first on the P wave components, then on the S and/or ST 
wave components, in either order. Subsequently, the method of the present invention is applied to remove tool 
and casing arrivals. In the preferred embodiment, however, the tool and casing arrivals are removed first, fol- 
lowed by the removal of the reflected wave components, converted wave components and/or noise, first on 

so the P wave components, then on the S and/or ST wave components, in either order. Other orders of operation, 

however, can be performed. 

At step 302, the arrival time of each P wave component is detected. Due to the nonhomogeneity of the 
formation, the P wave component of each waveform arrives at a different time, relative to each other. In the 
preferr d embodiment, a threshold detection scheme is appli d to the waveforms to detect th ir resp ctive 
55 arrival times. However, other arrival time detection schemes, e.g., a zero-crossing det ctor having an ampli- 
tude threshold requirement, can be employed. 

At step 303, the waveforms are preferably time shifted such that their P wave component arrival times 
are aligned. Alignment of the P wav component arrival times is shown with reference to Figure 4. Alternatively, 
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the P wave components arrival times can be detected and stored for use in conjunction with the remaining 
steps described hereinafter, thereby eliminating the need to align the wave components. 

At step 304, a time window is selected. Preferably, a 500 microsecond (us) time window is employed. Given 
the sampling rate of 10 us, the window will contain 50 sampled points for the P wave components of each wa- 
veform. . 

The time window preferably captures the entire P wave component. The P wave component is typically 
300 to 400 |is long. Additionally, the time window preferably commences 50 us before the detected P wave 
component arrival. This window position insures that the entire P wave component is captured for processing, 
as the P wave component's arrival time may be too small for detection and/or buried by noise. Thus, the time 
window of 500 us is preferably employed. Other time window values, either shorter or longer, can also be used. 

At step 305, the eigenvectors are obtained. The n-samples of each P wave component can be represented 
as a P wave component vector in n-space. The Karhunen-Loeve transformation generates eigenvectors or- 
thogonal to each other, derived from the statistics of the P wave component vectors. Eigenvectors are gener- 
ated based on a correlation matrix, CM, defined by the following equation: 



CM= lim f(lSjST) 

It— > oo k ^ j-l J J/ 



20 where 

k represents the number of vectors in the data set; 

Sj represents the vector, and 

Sj T represents the transpose of the vector. 
1 The correlation matrix definition is based on the statistical characteristics of the vectors given an infinite 
25 data set. As the data set of sonic logging is limited to m waveforms, the correlation matrix is approximated by 

the following equation: 



30 



CM=i(iSjSX) 



where 

m represents the number of waveforms in the data set; 
Sj represents the wave component vector; and 
35 S] T represents the transpose of the wave component vector. 

Given the P wave component vectors, the correlation matrix CM is calculated. The eigenvectors are ob- 
tained from the correlation matrix according to the following equation: 

CMOi = >LjOj 

where 

40 Oj represents the ith eigenvector; and 

Xi represents the ith scalar for the ith eigenvector. 

The scalar for the ith eigenvector is commonly referred to as its eigenvalue. The eigenvalues are preferably 
ranked according to magnitude, the first eigenvalue being the largest. As discussed above, the first eigenvalue 
represents the average of the wave component vectors in the first eigenvector's direction. The second through 

45 mth eigenvalues represent the average variation of the wave component vectors in the second through mth 
eigenvector's direction, respectively. 

In the preferred embodiment, the m eigenvectors are obtained using the IMSL software package, available 
from IMSL, Incorporated, Customer Relations, 2500 ParkWest Tower One, 2500 CityWest Boulevard, Houston, 
Texas, 77042-3020. However, other software packages are available and can be used for calculating the cor- 

so relation matrix, the eigenvectors and their corresponding eigenvalues. As will be appreciated by those skilled 
in the art, the eigenvectors contain the same number of points as there are samples in the P wave component. 

At step 306, the first eigenvector is correlated to each P wave component in the data set, thereby obtaining 
a correlation factor for each P wave component. In the preferred embodiment, the correlation factors are ob- 
tained by taking the inner product of the first eigenvector and each P wave component. The inner product is 

55 obtained by multiplying each first eigenvector point with each P wave component sample and cumulating the 
results to obtain a scalar value. Obtaining the correlation factor of the first eigenvector to the individual P wave 
components can be represented by the following equation: 
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a, = ]T E](i)S(i)j forj=ltom 
i = i 

5 

where 

aj represents the correlation factor for the jth waveform; 

E^i) represents the ith sample of the first eigenvector; 

S(i)j represents the ith sample of the P wave components contained in the time window; 

10 n represents the number of samples; and 

m represents the number of waveforms contained in the data set. 

At step 307, a correlated first eigenvector is removed from its respective P wave components, preferably 
by subtraction. In the preferred embodiment, the correlated first eigenvector is obtained for each of the m wa- 
veforms by multiplying each of the first eigenvector points by the correlation factor. Removing the correlated 

15 first eigenvector from the individual P wave component for each waveform can be represented by the following 
equation: 

S R (i)j = S(i)j - ajE^i) fori = 1 to n; j = 1 to m 

where 



S R (i)j represents the ith sample of the residual samples in the time window for the jth waveform; 

20 S(i)j represents the ith sample of the P wave component for the jth waveform; 

E^i) represents the ith sample of the first eigenvector; 

n represents the number of samples; and 

m represents the number of waveforms in the data set. 



Removing the correlated first eigenvector from the P wave component yields a residual for each of the 
25 waveforms. The residual comprises reflected wave components, converted wave components and noise, as 
well as a residual P wave component. 

As stated above, the first eigenvector captures about 86% of the wave component. Should greater accu- 
racy be desired, additional eigenvectors are employed, as set forth in steps 308 through 310. Otherwise, only 
steps 311 through 313 are followed. 
30 At step 308, the second eigenvector is correlated to each residual wave component, preferably as dis- 

cussed above with reference to step 306. Obtaining the correlation factor of the second eigenvector to each 
residual wave component can be represented by the following equation: 

n 

35 b J= X E 2 (i)S R (i)j forj = 1 torn 

i = i 

where 



40 bj represents the correlation factor for the jth waveform; 

E 2 (i) represents the ith sample of the second eigenvector; 

S R (i)j represents the ith sample of the residual wave components for the jth waveform; 

n represents the number of samples; and 

m represents the number of waveforms contained in the data set. 



At step 309, the correlated second eigenvector is removed from its respective residual wave components, 
preferably as discussed above with reference to step 307. Removing the correlated second eigenvector from 
the respective residual wave components can be represented by the following equation: 

S' R (i)j = S R (i)j - bjE 2 (i) for i = 1 to n; j = 1 to m 

where 



50 S' R (i)j represents the ith sample of the second residual for the jth waveform; 

S R (i)j represents the ith sample of the first residual for the jth waveform; 

E 2 (i) represents the ith sample of the second eigenvector; 

n represents the number of samples; and 

m represents the number of waveforms in the data set. 



As stated above, the first and second eigenvectors capture about 92% of the wave component. Should 
greater accuracy be desired, additional eigenvectors may employed at step 310. In the preferred embodiment, 
step 310 parallels steps 308 and 309, wherein the jth eigenv ctor is correlated to the 0-1 )th residual wave com- 
ponent sampl s, and the correlated jth eigenvector is subtracted from the most recent residual wave compo- 
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nent samples calculated. 

At step 311, the residual wave components are removed from the aligned P wave components, preferably 
by subtraction. In this case where only the first eigenvector was employed (to step 307), the P wave compo- 
nents are about 86% free of undesired wave components. In this case, removing the residual wave compo- 
nents from the aligned P wave components for each waveform can be represented by the following equation: 

S c (i)j = S(i)j - S R (i)j for i = 1 to n; j = 1 to m 

where 



S c (i)j represents the ith sample of the jth wave component substantially free of reflected wave compo- 
nents, converted wave components and/or noise; 
w S(i)j represents the ith sample of the jth wave component; 

S R (i)j represents the ith sample of the jth residual wave component; 
n represents the number of samples; and 

m represents the number of waveforms. 



In the case where the first and second eigenvectors were employed (to step 309), the P wave compon nts 
is are about 92% free of undesired wave components. In this case, removing the residual wave components from 
the aligned P wave components for each waveform can be represented by the following equation: 

S c (i)j = S(i)j - S' R (i)j for i = 1 to n; j = 1 to m 

where 



S c (i)j represents the ith sample of the jth wave component substantially free of reflected wave com- 

20 ponents, converted wave components and/or noise; 

S(i)j represents the ith sample of the jth wave component; 

S' R (i)| represents the ith sample of the jth residual wave component; 

n represents the number of samples; and 

m represents the number of waveforms. 



25 At step 312, the aligned waveforms are time shifted to return them to their original time positions, relative 

to each other. As shown with reference to Figure 5, the waveforms are shown with the first residual removed, 
as a result of step 307. As shown with reference to Figure 6, the waveforms are shown with the first and second 
residuals removed, as a result of step 309. 

As set forth in step 31 3, in order to remove undesired wave components from the shear (S) and/or Stoneley 
30 (ST) wave components, steps 302 through 312 are repeated. It is to be understood that "S" and "ST" would 
replace the "P" term through out Figures 3a and 3b. 

Turning now to Figure 7, a second embodiment of the method of the present invention for substantially 
removing reflected wave components, converted wave components and/or noise from acquired sonic data em- 
ploying an approximation to the Karhunen-Loeve transformation is described. In the second embodiment of 
35 the present invention, the first and second eigenvectors are approximated, thereby replacing steps 305 through 
311 of the first embodiment, as shown in Figures 3a and 3b. 

At step 701, a template of the P wave components is obtained. In the preferred embodiment, the template 
is obtained by stacking the P wave components. Stacking preferably involves cumulating corresponding wa- 
veform samples, located in the predetermined time window established at step 305, and dividing these cumu- 
40 lated samples by the number of waveforms cumulated. 

The template of step 701 represents the average value at each sampled point without reflected wave com- 
ponents, converted wave components and noise therein. As reflected wave components, converted wave com- 
ponents and noise are mathematically random, these components statistically average out to zero. P waves, 
on the other hand, are coherent, thereby allowing the stacking procedure to produce non-zero results where 
45 non-zero P wave components are present. 

It should be noted that the time window may, either by design or otherwise, contain a portion of the S wave 
component. Due to the fact that the S wave component arrival times vary according to formation type, the S 
wave components should not all line up when the P wave components are aligned (step 303). Thus, any S 
wave components captured by the time window should not be coherent and tend to average out to zero after 
so stacking. 

At step 702, the template is normalized. In the preferred embodiment, a normalization factor for the tem- 
plate is obtained by cumulating the square of each sample point in the template and then taking the squar 
root of the cumulation. The normalization factor can be represented by the following equation: 

55 
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i = i 

5 

where 

N represents the normalization factor; 

T(i) represents the ith sample of the template; and 

10 n represents the number of template samples. 

The template is normalized by dividing each template point by the normalization factor. The normalization 
of the template can be represented by the following equation: 

T N (i) = ^ for i = 1 to n 
N 

15 where 

T N (i) represents the ith sample of the normalized template; 
T(i) represents the ith sample of the template; 

N represents the normalization factor; and 

n represents the number of template samples. 

20 The template represents the average P wave component for the m waveforms in the data set. The nor- 

malized template represents an approximation of the first eigenvector. 

At step 703. the normalized template is correlated to each P wave component. In the preferred embodi- 
ment, a correlation factor for each P wave component is obtained by taking the inner product of the normalized 
template and the respective P wave component. The inner product is preferably obtained by multiplying each 

25 normalized template sample with the corresponding P wave component sample and cumulating the results to 
obtain a scalar value. Obtaining the correlation factor of the normalized template to the P wave components 
can be represented by the following equation: 



30 



a j = £ T N (i)S(i)j for j = 1 to m 

i=i. 



where 

35 a, represents the correlation factor for the jth waveform; 

Tfg(i) represents the ith sample of the normalized template; 
S(i) represents the ith sample of the jth P wave component; 

n represents the number of samples; and 

m represents the number of waveforms in the data set. 

40 At step 704, the correlated template is removed from its respective P wave components, preferably by sub- 

traction. In the preferred embodiment, the correlated template is obtained for each of the m waveforms by mul- 
tiplying the normalized template by the correlation factor for the respective waveform. Removing the correlated 
template from the respective P wave component can be represented by the following equation: 

S R (i)j = S(i)j - ajT N (i) for i = 1 to n; j = 1 to m 

45 where 

S R (i)j represents the ith sample of the jth residual; 

S(i)j represents the ith sample of the jth P wave component; 

a] represents the correlation factor for the jth wave component; 

T N (i) represents the ith sample of the normalized template; 
50 n represents the number of samples; and 

m represents the number of waveforms in the data set. 

Removing the normalized template from the respective P wave component yields a residual for each of the 
waveforms. The residual comprises reflected wave components, converted wave components and noise, as 
well as a residual P wave component. 
55 As with the first eigenvector, the approximation to the first eigenvector also captures about 86% of the 

wave component. Should greater accuracy be desired, an approximation to the second eigenvector can em- 
ployed, as set forth in steps 705 through 707. Otherwise, only step 708 is followed. 

At step 705, an approximation to the second igenvector is obtained. In th preferred embodiment, the 
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approximation to the second eigenvector is obtained by taking the Hilbert transform of the normalized template. 
The Hilbert transform is well known in the art and need not be discussed in detail herein. As will be appreciated 
in the art, the Hilbert transform of the normalized template yields a normalized Hilbert template. 

At step 706, the normalized Hilbert template is correlated to each residual wave component, preferably as 
5 discussed above with reference to step 703. Obtaining the correlation factor of the normalized Hilbert template 
to each residual P wave components can be represented by the following equation: 

n 

bj = X T H(i)SR(i)j for j = 1 to m 
10 t = l 



where 

bj represents the correlation factor for the jth waveform; 

15 T H (i) represents the ith sample of the Hilbert template; 

S R (i) represents the ith sample of the jth residual wave component; 

n represents the number of samples; and 

m represents the number of waveforms in the data set. 

At step 707, the correlated Hilbert template is removed from the residual wave components for which it 
20 was correlated, preferably as discussed above with reference to step 704. Removing the correlated Hilbert tem- 
plate from the respective residual wave components can be represented by the following equation: 

S' R (i)j = S R (i)j - bjT H (i) for i = 1 to n; j = 1 to m 

where 



S' R (i)j represents the ith sample of the jth second residual; 

25 S R (i)j . represents the ith sample of the jth first residual; 

bj represents the Hilbert correlation factor for the jth wave component; 

T H (i) represents the ith sample of the Hilbert template; 

n represents the number of samples; and 

m represents the number of waveforms in the data set. 



At step 708, the residual wave components are removed from the aligned P wave components, preferably 
by subtraction. In this case where only an. approximation to the first eigenvector was employed (to step 704), 
the P wave components are about 86% free of undesired wave components. In this case, removing the residual 
wave components from the aligned P wave components for each waveform can be represented by the following 
equation: 

S c (i)j = S(i)j - S R (i)j fori = 1 ton;j = 1 to m 

where 



S c (i)j represents the ith sample of the jth wave component substantially free of reflected wave compo- 
nents, converted wave components and/or noise; 
S(i)j represents the ith sample of the jth wave component; 
40 S R (i)j represents the ith sample of the jth first residual wave component; 
n represents the number of samples; and 

m represents the number of waveforms. 



In the case where approximations to the first and second eigenvectors were employed (to step 707), the 
P wave components are about 92% free of undesired wave components. In this case, removing the residual 
45 wave components from the aligned P wave components for each waveform can be represented by the following 
equation: 

S c (i)j = S(i)j - S' R (i)j for i = 1 to n; j = 1 to m 

where 



S c (i)j represents the ith sample of the jth wave component substantially free of reflected wave com- 

so ponents, converted wave components and/or noise; 

S(i)j represents the ith sample of the jth wave component; 

S' R (i)j represents the ith sample of the jth second residual wave component; 

n represents the number of samples; and 

m represents the number of waveforms. 



Thereafter, the aligned waveforms are time shifted to return them to their original time positions (step 31 2) 
and/or the steps are repeated to remove reflected wave components, converted wave components and/or 
noise from the S and/or ST wave components (step 313). 

In both the first and second embodiments, the time window for capturing the shear wave component is 
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preferably at least as wide as for the P wave component. The time window for capturing the Stoneley wave 
component, on the other hand, is preferably at least as wide as from the Stoneley arrival to the end of the re- 
cord d waveform. As is known in the art, the Stoneley wave component is typically longer than the time for 
which the waveform is recorded. In the preferred embodiment, after the tool and casing arrivals have b en 

5 filtered, it is preferable to filter the P wave component, the Stoneley wave component, and then the S wave 
component. The P is preferably first, due to its earliest arrival time, and the Stoneley is next, due to its dominant 
amplitude relative to the S wave component. 

As is known in the art, boreholes which are selected for production logging are typically lined with a casing 
and cemented in place in order to isolate zones in the formation. A good cement seal between casing and for- 

10 mation allows the casing and cement seal to be substantially transparent the tool's acoustic pulses. However, 
when the casing is poorly bonded to the formation, either due to an initially poor cementation job or through 
degradation over time, a portion of the tool's acoustic pulse is captured by the casing, resonating therein and 
detected by to the tool as casing arrivals. Casing arrivals, which typically appear before the P wave component, 
are recorded at substantially the same time over the course of the affected area. 

15 Tool arrivals are the product of the acoustic wave traveling from the transmitter, through the tool, to the 

receiver(s). As the spacing between the transmitter and receiver remains constant, tool arrivals are always 
recorded at substantially the same time, for a given transmitter-receiver pair, throughout the logging process. 

Turning now to Figure 8, a first embodiment of the present invention for substantially removing casing and 
tool arrivals is illustrated. After the data has been acquired at step 801 , the first eigenvector is obtained at step 

20 802 for the entire waveform, preferably as set forth at step 305. Due to the fact that the casing and tool arrivals 
appear at substantially the same time, it is not necessary to align the arrivals. Additionally, as the P, S and ST 
wave components will typically appear randomly relative to the casing and tool arrivals, it is not necessary to 
employ a time window. 

At step 803, the first eigenvector is correlated to each waveform in the data set, thereby obtaining a cor- 
25 relation .factor for each waveform, preferably as set forth at step 306. 

At step 804, a correlated first eigenvector is removed from its respective waveforms, preferably as set forth 
at step 307. Removing the correlated first eigenvector from each waveform yields a residual for each of the 
waveforms. The residual comprises reflected wave components, converted wave components and noise, as 
well as P, S and ST wave components. The residuals, however, do not contain an appreciable amount of casing 
30 or tool arrivals. However, should greater accuracy be desired, additional eigenvectors are employed, as set 

forth in steps 805 through 807. 

At step 805, the second eigenvector is correlated to each residual waveform, preferably as discussed 

above with reference to step 306. 

At step 806, the correlated second eigenvector is removed from its respective residual waveforms, pre- 

35 ferably as discussed above with reference to step 309. 

Additional eigenvectors may be employed at step 807, preferably as set forth at step 310. Thereafter, the 
residual waveform contains the P, S and ST wave components substantially free of casing and tool arrivals, 
and the reflected wave components, converted wave components and/or noise may be removed as explained 
above with reference to Figures 3 and/or 7. 

40 With reference to Figure 9, an acquired sonic data set is illustrated having tool arrivals 901 through 904 

which substantially mask the shear wave components. Employing the first eigenvector and the method dis- 
cussed above, the tool arrivals are removed, allowing the shear waves to be more readily illustrated, as shown 

with reference to Figure 10. 

Alternatively, in a second embodiment of the present invention for substantially removing casing and tool 
45 arrivals, approximations to the first and second eigenvectors may be employed. Thus, steps 802 through 807 
would be replaced with steps 701 through 707. 

In the preferred embodiment, two eigenvectors, or their approximations, are employed for removing casing 
arrivals, while one eigenvector, or its approximation, is preferably employed for removing tool arrivals. 

In another embodiment of the present invention, the reflection coefficient of a bed boundary in the forma- 
so tion can be estimated for each wave component. The reflection coefficient is estimated by dividing the energy 
content of the reflected wave component, obtained from the residual waveform, by the energy content of the 
transmitted wave component, obtained from the filtered wave component. In the preferred embodiment, the 
energy content of the reflected and converted wave components are calculated by cumulating the square the 
amplitude of each sample. Other m thods of calculating the en rgy content of the reflected and transmitted 
55 wave components are well known to thos skilled in the art and can be employ d in the alternative. 

The method of the present invention for estimating a reflection coefficient of a bed boundary in the for- 
mation repres nts an alternative to methods known in the art. The prior art calculated the reflection coefficient 
according to the following equation: 
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P2V2 ~ P1V1 

p 2 V 2 + PiV! 

where 

the transmitted wave travels through layer 1 and is reflected by layer 2; 
5 V | represents the velocity of the respective layers; and 
p represents the density of the respective layers. 

The velocity and density measurements of the layers were typically acquired in the prior art from a sonic 
and density log, respectively. Thus, both a sonic and a density logging tool have to be employed. The above 
method of the present invention for estimating a reflection coefficient based on transmitted and reflected wave 
10 components requires only the use of a sonic tool. Thus, substantial savings of both time and money will be 
realized by employing this method of the present invention. 

In yet another embodiment of the present invention, the lithology of the formation traversing the borehole 
can be identified. Acoustic pulses which travel through like media tend to exhibit like characteristics. I have 
found that plotting wave components based on their projection on the eigenvectors results in a cluster of wave 
15 components based on the media through which they traveled. As an example, Figure 11 shows a plurality of 
P wave components plotted in a two-dimensional plane. The abscissa (x-axis) is the projection of the P wave 
component on the first eigenvector, or its approximation, while the ordinate (y-axis) is the projection of the P 
wave component on the second eigenvector, or its approximation. 

It turns out that these points tend to cluster based on lithology type. For example, points obtained from 
20 wave components traveling through limestone tend to congregate in one area of the plane, while points ob- 
tained from wave components traveling through shaley sand tend to congregate in another area. Figure 11 
identifies four types of lithologies, e.g., shale, limestone and two different sand types, using the projection of 
the P wave component on the first and second eigenvectors for the acquired sonic data of Figure 2. 

The projection on the first eigenvector is the factor a J( calculated at step 306 or 703, correlating the P wave 
25 component to the first eigenvector or its approximation, respectively. The projection on the second eigenvector 
is the factor bj calculated at step 308 or 706, correlating the P wave component to the second eigenvector or 
its approximation, respectively. 

Given a plurality of wave components depicted in this manner, a database can be established to identify 
the lithology of the formation traversing the borehole given its coordinates in a plane whose axes are the com- 
30 ponent's projections on their eigenvectors. 

It is noted that, while general lithology can be identif ied using the projection of the respective P wave com- 
ponents on the first eigenvector, more specific lithology can be identified using the projection of the respective 
P wave component on the second eigenvector. It is further noted that the database could be any dimension, 
up to the n-dimensional space, based on the n samples of the wave component. 
35 |t is noted that any combination of wave components and/or eigenvectors may be employed in establishing 

the database. For example, the same wave component may be employed with its projections on differ nt ei- 
genvectors (as discussed above), different wave components may be employed with their respective projec- 
tions on the same eigenvector, or different wave components may be employed with their respective projec- 
tions on different eigenvectors. 
40 in still another embodiment of the present invention, data compression can be accomplished using the ei- 

genvectors as calculated above. Recalling that Karhunen-Loeve rotates the n-dimensional coordinate system 
to a new coordinate system, the wave component vectors can substantially be expressed in terms of their pri- 
mary eigenvectors. For example, if the 4th through 50th eigenvalues are substantially equal to zero, each wave 
component vector can substantially be described using only their projection on the first three eigenvectors. 
45 Assuming the wave component vectors can substantially be described using the f irst three eigenvectors, 

reconstruction factors A 1f A 2 and A 3 are defined as follows: 

Aj = EjSj 
Bj = E; Sj 

50 J 3 J 

where 

Aj represents a scalar reconstruction factor for the 1st eigenvector, jth waveform; 
Bj represents a scalar reconstruction factor for the 2nd eigenvector, jth waveform; 
Cj represents a scalar reconstruction factor for the 3rd eigenvector, jth waveform; 
55 El represents the transpose of the ith eigenvector; and 

Sj represents the jth wave component vector. 

It is noted that A jt Bj and Cj are the vectors obtained from correlating the first, second and third eigenvector 
to the individual wave components, respectively, i.e., the vector obtained from the individual aj's, b/s and c/s 
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calculated as set forth at steps 306,308 and 310, respectively. 

The three scalar reconstruction factors are retained for each wave component vector, as well as the three 
eigenvectors. For illustrative purposes, assume each wave component vector has 50 samples and the data 
set has 10,000 waveforms. Without compression, the data set would require 1500k values for storage (50 sam- 
ples per wave component * 3 wave components per waveform * 10k waveforms). By the data compression 
method herein described, only 90.45k values are required (3 wave components per waveform * 3 reconstruc- 
tion factors per wave component * 10k waveforms + 3 wave components per waveform * 3 eigenvectors per 
wave component * 50 samples per eigenvector). The compression ratio of the present invention is thus about 
1 7 to 1 . 

In order to reconstruct the original wave components, the following equation is preferably employed: 

Sj = AjEi + BjE 2 + CjE 3 forj = 1 to m 

The above techniques have been described in detail with reference to sonic data acquired with one trans- 
mitter and one receiver. These techniques are also applicable to sonic data acquired by an array sonic tool, 
having at least one transmitter and a plurality of receivers, as well as a borehole compensation tool, having 
at least two transmitters and at least two receivers. Such tool configurations are well known and need not be 
described in detail herein. 

Although illustrative embodiments of the present invention have been described in detail with reference 
to the accompanying drawings, it is to be understood that the invention is not limited to those precise embodi- 
ments. Various changes or modifications may be effected therein by one skilled in the art without departing 
from the scope or spirit of the invention. 



Claims 

1. Arpethod of filtering acquired sonic well logging data comprising a data set of m waveforms each including 
a component comprising n samples characterised as a vector, the method comprising obtaining an eigen- 
vector based on said vectors and using the eigenvector to filter out unwanted components contained in 
the data, characterised in that the data comprise sonic data, each component comprises a formation wave 
component digitized into said n samples and said vectors comprise formation wave component vectors, 
the method further comprising the steps of: 

obtaining a first eigenvector based on said formation wave component vectors; correlating said first 
eigenvector to each of said wave component vectors, thereby obtaining a correlated first eigenvector for 
each of said wave component vectors; 

removing said correlated first eigenvectors from their respective wave component vectors, thereby 
obtaining a first residual component for each of said wave component vectors; and 

removing said first residual components from their respective formation wave components, there- 
by substantially removing the reflected wave components, converted wave components and/or noise from 
said formation wave component. 

2. The method of claim 1 , wherein the step of removing said first residual components from their respective 
formation wave components is based on the following equation: 

S c (i)j = S(i)j - S R (i)j fori = 1 to n; j = 1 torn 

where 

S c (i)j represents the ith sample of a formation wave component for the jth waveform whose reflected 
wave components converted wave components and/or noise has been substantially removed therefrom; 
S(i)j represents the ith sample of the formation wave component vector for the jth waveform; and 
S R (i)j represents the ith sample of the first residual component for the jth waveform. 

3. The method of claim 1 or 2, wherein said formation wave component is selected from a compressional 
wave component, a shear wave component or a Stoneley wave component. 

4. The method of any preceding claim, 

said method further including the steps of: 

calculating the energy cont nt of at least a portion of said reflected wave component removed from 

said formation wave component; 

calculating the energy content of at least a portion of said formation wave component; and 
dividing said calculated energy content of said reflected wave component by said calculated en rgy 

content of said formation wave component, thereby obtaining an estimate of a reflection coeff ici nt of a 
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bed boundary in the formation. 

A method of filtering acquired sonic well logging data comprising a data set of m waveforms each including 
a component comprising n samples characterised as a vector, the method comprising obtaining an eigen- 
vector based on said vectors and using the eigenvector to filter out unwanted components contained in. 
the data, characterised in that the data comprise sonic data, each component comprises a formation wave 
component digitized into said n samples and said vectors comprise formation wave component vectors, 
the method further comprising the steps of: 

obtaining a first eigenvector based on said formation waveform vectors; 

correlating said first eigenvector to each of said formation waveform vectors, thereby obtaining a 
correlated first eigenvector for each of said waveform vectors; 

removing said correlated first eigenvectors from their respective formation waveform vectors, 
thereby substantially removing the casing and/or tool arrivals from said formation waveforms. 

The method of any preceding claim, wherein said step of obtaining a first eigenvector includes the steps 
of: 

transposing each of said formation wave component vectors; 

obtaining a correlation matrix based on said formation wave component vectors and said transpose 
thereof; and 

obtaining a first eigenvector based on said correlation matrix. 

The method of any preceding claim, wherein the step of removing said correlated first eigenvector from 
its respective wave component vector is based on the following equation: 

S R (i)i = S(i)j - ajE^i) for i = 1 to n; j = 1 to m 

where 

S R (i)j represents the ith sample of the first residual component for the jth waveform; 
S(i)j represents the ith sample of the formation wave component vector for the jth waveform; 
aj represents a correlation factor correlating said first eigenvector to the jth formation wave com- 
ponent; and 

E^i) represents the ith sample of said first eigenvector. 

A method of filtering acquired sonic well logging data, the sonic data comprising a data set of m waveforms, 
each of the waveforms including a formation wave component digitized into n samples, said filtering meth- 
od to substantially remove reflected wave components, converted wave components and/or noise from 
the formation wave component, said method comprising the steps of: 

stacking the digitized formation wave components, thereby obtaining a template of the formation 
wave components; 

normalizing said template, thereby obtaining an average value of said formation wave components; 

correlating said normalized template to each of said formation wave components, thereby obtaining 
a correlated template for each of said formation wave components, said template having n samples; 

removing said correlated template from the respective formation wave components, thereby ob- 
taining a first residual component for each of said formation wave components; and 

removing said first residual components from their respective formation wave components, there- 
by substantially removing the reflected wave components, converted wave components and/or noise from 
said formation wave component. 

The method of claim 8, wherein the step of removing said first residual components from their respective 
formation wave components is based on the following equation: 

S c (i)j = S(i)j - S R (i)j for i = 1 to n; j = 1 torn 

where 

S c (i)] represents the ith sample of a formation wave component for the jth waveform whose reflected 
wave components, converted wave components and/or noise has been substantially removed therefrom; 
S(i)j represents the ith sample of the formation wave component for the jth waveform; and 
S R (i)j represents the ith sample of the f irst residual component for the jth waveform. 

0. The method of claim 8 or 9, wherein said formation wave component is selected from a compressional 
wave component, a shear wave component or a Stoneley wave component. 

1. The method of any of claims 8 to 10, said method further including the steps of: 
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calculating the energy content of at least a portion of said reflected wave component removed from 

said formation wave component; 

calculating the energy content of at least a portion of said formation wave component; and 
dividing said calculated energy content of said reflected wave component by said calculated energy 

content of said formation wave component, thereby obtaining an estimate of a reflection coefficient of a 

bed boundary in the formation. 

12. A method of filtering acquired sonic well logging data, the sonic data comprising a data set of m formation 
waveforms each of the waveforms being digitized into n samples, said method comprising the steps of: 

stacking the formation waveforms, thereby obtaining a template of the formation waveforms; 

normalizina said template. Thereby obtaining an average value of said formation waveforms; 

correlating said normalized template to each of said formation waveforms, thereby obtaining a cor- 
related tempiate for each of said formation waveforms, said template having n samples; and 

removing said correlated template from the respective formation waveforms, thereby substantially 
removing the casing and/or tool arrivals from said formation waveforms. 

13. The method of any of claims 8-12, said step of normalizing said template comprising the steps of: 

cumulating the square of each template sample; 

obtaining the square root of said cumulation, thereby obtaining a normalization factor; and 
dividing each of said template samples by said normalization factor, thereby normalizing said tem- 
plate. 

14. The method of any of claims 8-13, wherein the step of removing said correlated template from the respec- 
tive formation waveforms is based on the following equation: 

S c (i)j = S(i)j - ajT N (i) for i = 1 to n; j = 1 to m 

where 

S c (i)j represents the ith sample of the formation waveform for the jth waveform whose casing and/or 
tool arrival has been substantially removed therefrom; 

S(i)j represents the ith sample of the formation waveform for the jth waveform; 

aj represents the correlation factor correlating the normalized template to the jth formation wave- 
form; and 

T N (i) represents the ith sample of the normalized template. 



Patentanspruche 

1 . Ein Verfahren zum Filtern gewonnener akustischer Bohrloch-Logdaten, umfassend einen Datensatz von 
m-Wel!enformen mit jeweils einer Komponente, die n als ein Vektor charakterisierte Abtastwerte umfafit, 
welches Verfahren das Gewinnen eines Eigenvektors umfaftt, basierend auf den Vektoren und unter Ver- 
wendung des Eigenvektors zum Ausfiltern ungewunschter in den Daten enthaltener Komponenten, da- 
durch gekennzeichnet, da& die Daten akustische Daten umfassen, jede Komponente eine Formations- 
wellenkomponente umfalit, die in die n-Abtastwerte digitalisiert ist und die Vektoren Formationswellen- 
komponentenvektoren umfassen, welches Verfahren ferner die Schritte umfa&t: 

Gewinnen eines ersten Eigenvektors, basierend auf den Formationswellenkomponentenvektoren; 

Korrelieren des ersten Eigenvektors mit jedem der Wellenkomponentenvektoren, wodurch ein kor- 
relierter erster Eigenvektor fur jeden der Wellenkomponentenvektoren gewonnen wird; 

Entfernen der korrelierten ersten Eigenvektoren aus ihren jeweiligen Wellenkomponentenvekto- 
ren, wodurch eine erste Restkomponente fur jeden der Wellenkomponentenvektoren gewonnen wird; und 
Entfernen der ersten Restkomponente aus den jeweiligen Formationswellenkomponenten, wodurch re- 
flektierte Wellenkomponenten, konvertierte Wellenkomponenten und/oder Rauschen aus der Formati- 
onswellenkomponente im wesentlichen entfernt werden. 

2. Das Verfahren nach Anspruch 1, bei dem der Schritt des Entfernens der ersten Restkomponenten aus 
ihren jeweiligen Formationswellenkomponenten auf der folgenden Gleichung basiert: 

S c (i)j = S(i)j - S R (i)j fur i = 1 bis n;j = 1 bism 

worin 

S c (i)jdeni-tenAbtastw rt einer Formationswellenkomponente fur die j-te Wellenform reprasentiert, 
deren reflektierte Wellenkomponenten, konvertierte Wellenkomponenten und/oder Rauschen aus ihr im 
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wesentlichen entfernt ist; 

S(i)j den i-ten Abtastwert des Formationswellenkomponentenvektors fur die j-te Wellenform repra- 

sentiert; und 

S R (i)j den i-ten Abtastwert der ersten Restkomponente fur die j-te Wellenform reprasentiert. 

Das Verfahren nach Anspruch 1 Oder 2, bei dem die Formations-wellenkomponente ausgewahlt wird aus 
einer Kompressionswellenkomponente, einer Scherwellenkomponente oder einer Stoneley- 
Wellenkomponente. 

Das Verfahren nach einem der vorangehenden Anspruche, welches Verfahren ferner die Schritte umfafit: 
Berechnen des Energieinhalts von zumindest einem Abschnitt der ref lektierten Welle nkom pone n- 

te, entfernt aus der Formationswellenkomponente; 

Berechnen des Energieinhalts von zumindest einem Abschnitt der Formationswellenkomponente; 

und 

Dividieren des berechneten Energieinhalts der reflektierten Wellenkomponente durch den berech- 
neten Energieinhalt der Formationswellenkomponente, wodurch man eine Abschatzung eines Reflexions- 
koeff izienten einer Schichtgrenze in der Formation erhalt. 

Ein Verfahren zum Filtern gewonnener akustischer Bohrloch-Logdaten, umfassend einen Datensatz von 
m-Wellenformen, welche jeweils eine Komponente einschlie&en, umfassend n als ein Vektor charakteri- 
sierte Abtastwerte, welches Verfahren das Gewinnen eines Eigenvektors, basierend auf den Vektoren, 
und Verwenden des Eigenvektors zum Ausf iltern unerwunschter in den Daten enthaltener Komponenten 
umfafit, dadurch gekennzeichnet, dad die Daten akustische Daten umfassen, jede Komponente eine For- 
mationswellenkomponente umfa&t, die in die n-Abtastwerte digitalisiert ist und die Vektoren Formations- 
wellenkomponentenvektoren umfassen, welches Verfahren ferner die Schritte umfafit: 

Gewinnen eines ersten Eigenvektors, basierend auf den Formationswellenformvektoren; 

Korrelieren des ersten Eigenvektors mit jedem der Formationswellenformvektoren, wodurch man 
einen korrelierten ersten Eigenvektor fur jedes der Wellenformvektoren erhalt; 

Entfernen des korrelierten ersten Eigenvektors aus den jeweiligen Formationswellenformvektoren, 
wodurch die Auskleidungs- und/oder Sondeneinlaufe von den Formationswellenformen im wesentlichen 
entfernt werden. 

Das Verfahren nach einem der vorangehenden Anspruche, bei dem der Schritt des Gewinnens eines er- 
sten Eigenvektors die Schritte umfafit: 

Transponieren jedes der Formationswellenkomponentenvektoren; 

Gewinnen einer Korrelationsmatrix, basierend auf den Formationswellenkomponentenvektoren 

und deren Transposition; und 

Gewinnen eines ersten Eigenvektors, basierend auf der Korrelationsmatrix. 

Das Verfahren nach einem der vorangehenden Anspruche, bei dem der Schritt des Entfernens des kor- 
relierten ersten Eigenvektors aus seinem entsprechenden Wellenkomponentenvektor auf der folgenden 
Gleichung basiert: 

S R (i)j = S(i)j - ajE^i) f ur i = 1 bis n; j = 1 bis m 

worin 

S R (i)j den i-ten Abtastwert der ersten Restkomponente fur die j-te Wellenform reprasentiert; 
S(i)j den i-ten Abtastwert des Formationswellenkomponentenvektors fur die j-te Wellenform repra- 
sentiert; 

a } einen Korrelationsfaktor reprasentiert, der den ersten Eigenvektor mit der j-ten Formationswel- 
lenkomponente korreliert; und 

E^i) den i-ten Abtastwert des ersten Eigenvektors reprasentiert. 

Ein Verfahren zum Filtern gewonnener akustischer Bohrloch-Logdaten, welche akustischen Daten einen 
Datensatz von m- Wellenform en umfassen, wobei jede der Weltenformen eine in n-Abtastwerte digitali- 
sierte Formationswellenkomponente umfafit, welches Filterverfahren dazu dient, reflektierte Wellenkom- 
ponenten, konvertierte Wellenkomponenten und/oder Rauschen aus der Formationswellenkomponente 
im wesentlichen zu entfernen, welches Verfahren die Schritte umfafit: 

Stapeln der digital isierten Formationswellenkomponenten, wodurch eine Schablone der Formati- 
onswellenkomponenten gewonnen wird; 

Normalisieren der Schablone, wodurch ein Mittelwert der Formationswellenkomponenten gewon- 
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nen wird; 

Korrelieren der normalisierten Schablone mit jeder der Formationswellenkomponenten, wodurch 
eine korrelierte Schablone fur jede der Formationswellenkomponenten gewonnen wird, welche Schablone 
n-Abtastwerte besitzt; 

Entfernen der korrelierten Schablone aus den entsprechenden Formationswellenkomponenten, 
wodurch eine erste Restkomponente fur jede der Formationswellenkomponenten gewonnen wird; und 

Entfernen der ersten Restkomponenten aus ihren entsprechenden Formationswellenkomponen- 
ten, wodurch die reflektierten Wellenkomponenten, konvertierten Wellenkomponenten und/oder Rau- 
schen aus der Formationswellenkomponente im wesentlichen entfernt werden. 

Das Verfahren nach Anspruch 8, bei dem derSchrittdes Entfernen der ersten Restkomponente aus ihren 
entsprechenden Formationswellenkomponenten auf der folgenden Gleichung basiert: 

S c (i)j = S(i)j - S R (i)j fur i = 1 bisn;j = 1 bis m 

worin 

S c (i)j den i-ten Abtast wert einer Formationswellenkomponente f ur die j-te Wellenform reprasentiert, 
deren reflektierte Wellenkomponenten, konvertierten Wellenkomponenten und/oder Rauschen aus ihr im 
wesentlichen entfernt worden ist; 

Sfljj den i-ten Abtastwert der Formationswellenkomponente fur die j-te Wellenform reprasentiert; 

und 

S R (i) i den i-ten Abtastwert der ersten Restkomponente fur die j-te Wellenform reprasentiert. 

0. Das Verfahren nach Anspruch 8 Oder 9, bei dem die Formationswellenkomponente ausgewahlt wird aus 
einer Kompressionswellenkomponente, einer Scherwellenkomponente Oder einer Stonel y- 
Wellenkomponente. 

1. Das Verfahren nach einem der Anspruche 8 bis 10, welches Verfahren ferner die Schritte umfafit: 

Berechnen des Energieinhalts von zumindest einem Abschnitt der reflektierten Wellenkomponen- 
te, entfernt aus der Formationswellenkomponente; 

Berechnen des Energieinhalts von zumindest einem Abschnitt der Formationswellenkomponente; 

und 

Dividieren des berechneten Energieinhalts der reflektierten Wellenkomponente durch den berech- 
neten Energieinhalt der Formationswellenkomponente, wodurch eine Abschatzung des Reflexionskoef- 
fizienten einer Schichtgrenze in der Formation gewonnen wird. 

2. Ein Verfahren zum Filtern gewonnener akustischer Bohrloch-Logdaten, welche akustischen Daten einen 
Datensatz von m-Formationswellenformen umfafit, wobei jede Wellenform in n-Abtastwerte digitalisiert 
ist, welches Verfahren die Schritte umfaftt: 

Stapeln der Formationswellenformen, wodurch eine Schablone der Formationswellenformen ge- 
wonnen wird; 

Normalisieren der Schablone, wodurch ein Mittel wert der Formationswellenformen gewonnen wird; 

Korrelieren der normalisierten Schablone mit jeder der Formationswellenformen, wodurch eine kor- 
relierte Schablone fur jede der Formationswellenformen gewonnen wird, welche Schablone n-Abtastwer- 
te besitzt; und 

Entfernen der korrelierten Schablone aus den jeweiligen Formationswellenformen, wodurch 
Auskleidungs- und/oder Sondeneinlaufe von den Formationswellenformen im wesentlichen entfernt wer- 
den. 

3. Das Verfahren nach einem der Anspruche 8 bis 12, wobei der Schritt der Normalisierung der Schablone 
die Schritte umfafit: 

Kumulieren des Quadrats jedes Schablonenabtastwertes; Gewinnen der Quadratwurzel der Ku- 
mulation, wodurch ein Normalisationsfaktor gewonnen wird; und 

Dividieren jedes Schablonenabtastwertes durch den Normalisationsfaktor, wodurch die Schablone 
normalisiert wird. 

4. Das Verfahr n nach einem der Anspruch 8 bis 13, bei dem der Schritt des Entfernens der korreli rten 
Schablone aus den j weilig n Formationswellenformen auf der folgenden Gleichung basiert: 

S c (i) j = S(i)j - ajT N (i) fur i = 1 bis n; j = 1 bis m 

worin 

S c (i)j den i-ten Abtastwert der Formationswellenform fur die j-te Wellenform reprasentiert, dessen 
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Auskleidungs- und/oder Sondeneinlauf davon im wesentlichen entfernt worden ist; 

S(i)j den i-ten Abtastwert der Formationswellenform fur die j-te Wellenform reprasentiert; 

^ den Korrelationsfaktor reprasentiert, der die normalisierte Schablone mit der j-ten Formations- 
wellenform korreliert; und 

T N (i) den i-ten Abtastwert der normalisierten Schablone reprasentiert. 

Revendications 

1. Procede de filtrage des donnees de diagraphie acoustique acquises comprenant un ensemble de don- 
nees de m formes d'onde, chacune incluant une composante comprenant n echantillons caracterises 
comme vecteur, le procede comprenant I'obtention d'un vecteur propre base sur lesdits vecteurs et utili- 
sant le vecteur propre pour filtrer les composantes superflues contenues dans les donnees, caracterise 
en ce que les donnees comprennent des donnees acoustiques, chaque composante comprend une 
composante d'onde de formation numerisee en lesdits n echantillons et lesdits vecteurs comprennent des 
vecteurs de composantes d'onde de formation, le procede comprenant en outre les etapes de : 

obtention d'un premier vecteur propre base sur lesdits vecteurs de composantes d'onde de forma- 
tion ; 

correlation dudit premier vecteur propre avec chacun desdits vecteurs de composantes d'onde, 
d'ou I'obtention d'un premier vecteur propre correle pour chacun desdits vecteurs de composantes 
d'onde ; 

suppression desdits premiers vecteurs propres correles de leurs vecteurs de composantes d'onde 
respectifs, d'ou I'obtention d'une premiere composante residuelle pour chacun desdits vecteurs de 
composantes d'onde ; et 

suppression desdites premieres composantes residuelles de leurs composantes d'onde de forma- 
tion respectives, d'ou une suppression substantielle des composantes d'onde reflechie, des composantes 
d'onde convertie et/ou du bruit de ladite composante d'onde de formation. 

2. Procede selon ia revendication 1, dans lequel I'etape de suppression desdites premieres composantes 
residuelles de leurs composantes d'onde de formation respectives est base sur I'equation suivante : 

S c (i)j = S(i)j - S R (i) } pour i = 1 a n ; j = 1 a m 

ou 

S c (i)j represente le i* me echantillon d'une composante d'onde de formation pour la j* me forme d'onde 
dont les composantes d'onde reflechie, les composantes d'onde convertie et/ou le bruit en ont ete subs- 
tantieilement supprimes ; 

S(i)j represente le i* me echantillon du vecteur de composante d'onde de formation pour la j* me forme 

d'onde ; et 

S R (i)j represente le i* me echantillon de la premiere composante residuelle pour la j« me forme d'onde. 

3. Procede selon la revendication 1 ou la revendication 2, dans lequel ladite composante d'onde de formation 
est choisie a partir d'une composante d'onde de compression, d'une composante d'onde de cisaillement 
ou d'une composante d'onde de Stoneley. 

4. Procede selon I'une quelconque des revendications precedentes, ladite methode incluant en outre les eta- 
pes de : 

calcul du contenu energetique d'au moins une partie de ladite composante d'onde reflechie sup- 
primee de ladite composante d'onde de formation ; 

calcul du contenu energetique d'au moins une partie de ladite composante d'onde de formation ; 

et 

division dudit contenu energetique calcule de ladite composante d'onde reflechie par ledit contenu 
energetique calcule de ladite composante d'onde de formation, d'ou I'obtention d'une estimation d'un 
coefficient de reflexion d'un bord de strate dans la formation. 

5. Procede de filtrage des donnees de diagraphie acoustique acquises comprenant un ensemble de don- 
nees de m formes d'ondes, chacune incluant une composante comprenant n echantillons caracterises 
comme vecteur, le procede comprenant I'obtention d'un vecteur propre base sur lesdits vecteurs et utili- 
sant le vecteur propre pour filtrer les composantes superflues contenues dans les donnees, caracterise 
en ce que les donnees comprennent des donnees acoustiques, chaque composante comprend une 
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composante d'onde deformation numerisee en lesdits n echantillons et lesdits vecteurs comprennent des 
vecteurs de composantes d'onde de formation, le procede comprenant en outre les etapes de : 

obtention d'un premier vecteur propre base sur lesdits vecteurs de forme d'onde de formation ; 

correlation dudit premier vecteur propre a chacun desdits vecteurs de forme d'onde de formation, 
d'ou I'obtention d'un premier vecteur propre correle pour chacun desdits vecteurs de forme d'onde ; 

suppression desdits premiers vecteurs propres correles de leurs vecteurs de forme d'onde de for- 
mation respectifs, d'ou une suppression substantielle des arrivages de tubage et/ou d'outils desdites for- 
mes d'onde de formation. 

Procede selon I'une quelconque des revendications precedentes, dans lequel ladite etape d'obtention 
d'un premier vecteur propre inclut les etapes de : 

transposition de chacun desdits vecteurs de composantes d'onde de formation ; 

obtention d'une matrice de correlation basee sur lesdits vecteurs de composantes d'onde de for- 
mation et de ladite transposition de ceux-ci ; et 

obtention d'un premier vecteur propre base sur ladite matrice de correlation. 

Procede selon I'une quelconque des revendications precedentes, dans lequel I'etape de suppression dudit 
premier vecteur propre correle de son vecteur de composantes d'onde respectif est base sur I'equation 
suivante : 

S R (i)j = S(i)j - ajE^t) pour i = 1 a n ; j = 1 a m 

ou ( 
S R (i)j represente le i« me echantilion de la premiere composante residuelle pour la forme d onde ; 
S(i)j represente le echantilion de la composante d'onde de formation pour la j*™ forme d'onde ; 
aj represente un facteur de correlation mettant en correlation ledit premier vecteur propre avec la 
jeme composante d'onde de formation ; et 

E«i(i) represente le i* me echantilion dudit premier vecteur propre. 

Procede de filtrage des donnees de diagraphie acoustique acquises, les donnees acoustiques compre- 
nant un ensemble de donnees de m formes d'onde, chacune des formes d'onde incluant une composante 
d'onde de formation numerisee en n echantillons, ledit procede de filtrage pour supprimer substantielle- 
ment les composantes d'onde reflechie, les composantes d'onde convertie et/ou le bruit de la composante 
d'onde de formation, ledit procede comprenant les etapes de : 

empilement des composantes d'onde de formation numerisees, d'ou I'obtention d'un gabarit des 

composantes d'onde de formation ; 

normalisation dudit gabarit, d'ou I'obtention d'une valeur moyenne desdites composantes d'onde 

de formation ; t % 

correlation dudit gabarit normalise avec chacune desdites composantes d'onde de formation, d'ou 
I'obtention d'un gabarit correle pour chacune desdites composantes d'onde de formation, ledit gabarit 
ayant n echantillons ; 

suppression dudit gabarit correlee des composantes d'onde de formation, d'ou I'obtention d'une 
premiere composante residuelle pour chacune desdites composantes d'onde de formation ; et 

suppression desdites premieres composantes residuelles de leurs composantes d'onde de forma- 
tion respectives, d'ou une suppression substantielle des composantes d'onde reflechie, de composantes 
d'onde convertie et/ou du bruit de ladite composante d'onde de formation. 

Procede selon la revendication 8, dans lequel I'etape de suppression desdites premieres composantes 
residuelles de leurs composantes d'onde de formation respectives est basee sur I'equation suivante : 

S c (i)j = S(i)j - S R (i)j pour i = 1 a n ; j = 1 a m 

ou 

S c (i)j represente le i*"* echantilion d'une composante d'onde de formation pour la j e ™ forme d'onde 
dont les composantes d'onde reflechie, les composantes d'onde convertie et/ou le bruit en ont ete subs- 

tantiellement supprimes ; 

S(i)j represente le i* me echantilion de la composante d'onde de formation pour la forme d'onde ; 

S R (i)j repres nte le echantilion de la pr miere composante residuelle pour la j* me forme d'onde. 

). Procede selon la revendication 8 ou la revendication 9, dans lequel ladite composante d'onde de formation 
est choisie a partir d'une composante d'onde de compression, d'une composante d'onde de cisaillement 
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ou cTune composante d'onde de Stoneley. 

Procede selon Tune des revendications 8 a 10, ledit procede incluant en outre les etapes de : 

calcul du contenu energetique d'au moins une partie de ladite composante d'onde reflechie sup- 

primee de ladite composante d'onde de formation ; 

calcul du contenu energetique d'au moins une partie de ladite composante d'onde de formation ; 

et 

division dudit contenu energetique calcule de ladite composante d'onde reflechie par iedit contenu 
energetique calcule de ladite composante d'onde de formation, d'ou I'obtention d'une estimation d'un 
coefficient de reflexion d'une limite de strate dans la formation. 

Procede de f iltrage des donnees de diagraphie acoustique acquises, les donnees acoustiques compre- 
nant un ensemble de donnees de m formes d'onde de formation, chacune des formes d'onde etant nu- 
merisee en n echantillons, ledit procede comprenant les etapes de : 

empilement des formes d'onde de formation, d'ou I'obtention d'un gabarit des formes d'onde de 

formation ; 

normalisation dudit gabarit, d'ou I'obtention d'une valeur moyenne desdites formes d'onde de for- 
mation ; 

correlation dudit gabarit normalise pour chacune des formes d'onde de formation, d'ou I'obtention 
d'un gabarit correle pour chacune des formes d'onde de formation, ledit gabarit ayant n echantillons ; et 

suppression dudit gabarit correle des formes d'onde de formation respectives, d'ou une suppres- 
sion substantielle des arrivages de tubage et/ou d'outils desdites formes d'onde de formation. 

Procede selon I'une des revendications 8-12, ladite etape de normalisation dudit gabarit comprenant les 
etapes de : 

cumul du carre de chaque echantillon de gabarit ; 

obtention de la racine carr6e dudit cumul, d'ou I'obtention d'un facteur de normalisation ; et 
division de chacun desdits echantillons de gabarit par ledit facteur de normalisation, d'ou la nor- 
malisation dudit gabarit. 

Procede selon I'une des revendications 8-1 3, dans lequel I'etape de suppression dudit gabarit correle des 
formes d'onde de formation respectives est basee sur I'equation suivante : 

S c (i)j = S(i)j - ajT N (i) pour i = 1 a n ; j = 1 a m 

ou 

S c (i)j represente le i« me echantillon de la forme d'onde de formation pour la j* me forme d'onde dont 
I'arrivage de tubage et/ou d'outils en a ete substantiellement supprime ; 

S(i)j represente le echantillon de la forme d'onde de formation pour la j* me forme d'onde ; 

aj represente le facteur de correlation mettant en correlation le gabarit normalise avec la j* me forme 
d'onde de formation ; et 

T N (i) represente le i* me Echantillon du gabarit normalise. 
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